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ABSTRACT 


A parametric study of a body of revolution is conducted utilizing existing semi- 
empirical methods for the calculation of hydrodynamic coefficients. The geometry of 
the body is analyzed in non-dimensional length and volume parameters. The effects of 
varying the nose, mid-body, and base fractions of the body on the hydrodynamic 
coefficients are generated and illustrated graphically. Equations for the hydrodynamic 
coefficients are then determined from the non-dimensional parameters. The results 
can be used to evaluate fundamental maneuvering characteristics early in the design 


phase. 
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I. INTRODUCTION 


A. GENERAL CONCEPTS 


Ships are designed to operate in a hostile environment. Even 
during peacetime the ship must navigate the water and arrive at it's 
destination intact and ontime. The motion of the ship is affected by 
the presence of other ships, hazards to navigation, established traffic 
patterns, and the sea. The motion of the ship is directed and the 
track of the ship is plotted on the chart. The response of the ship to 
changes in speed and rudder commands is studied under standard 
assumptions of marine vehicle dynamics. The forces and moments 
subjected to the ship due to its relative motion in the water are the 
basis for this study. Figure 1 illustrates the basic concepts and 


definitions. 
Y Force | Vv 


N Moment r 






Lift Forces uo 


Drag Forces Velocities 


Rudder Angle of Attack 
Angle 


Figure 1. Overview of forces and velocities. [1] 











The motion of the ship can be described in the global sense and 
in the body fixed reference frame. The body fixed reference frame is 
utilized in this study. The ship's velocity along its x axis is u. The 
sideslip velocity is v and the angular velocity is r. 

When the water impacts the ship at an angle of attack, lift and 
drag forces are developed. These impart a lateral Y force and a 
pitching N moment. We can determine the lateral Y force resulting 
from a sideslip velocity v and an angular velocity r. We can 
determine the pitching N moment resulting from a slideslip velocity 
v and angular velocity r. In review of the effects of acceleration, 

v and; , these effects can be computed. 

A measure of the ship's stability in design is its ability to 
maintain or regain straight line motion when forces act on its shape. 
The linearized models for equations of motion utilize hydrodynamic 
derivities of the Y force and the N moment to describe their 
behavior. There are a number of mathematical methods for 
determining the hydrodynamic derivatives and LT Wolkerstofer [2] 
Summarized them in his thesis. In our work, the semi-empirical 
methods that utilize the geometric considerations of a body of 
revolution, typical of a modern submarine are best suited for our 


purpose. 


B. PROGRAM APPLICATION 


The body of revolution that we shall consider is a basic 


submarine shape. The nose is elliptical, the mid-body is cylindrical, 








and the base is conical. This is an approximate SUBOFF model [2] and 


it is also the shape of the slice pod which is illustrated in Figure 2. 


Figure 2. A body of revolution. 


This axisymmetric shape is a simple model but the concepts 
applied to this model can be broadened to fit any shape. Our goal is 
to take a shape that has been tested and to modify its parameters to 
see how geometric changes effect the hydrodynamic derivatives 


which affect the stability and turning characteristics of the platform. 








fi. DETERMINATION OF HYDRODYNAMIC COEFFICIENTS 


A. THE GENERAL PHYSICAL DESCRIPTION 


The body of revolution consists of three main sections. The 
forward section is an ellipsoid, the mid-body section is a cylinder and 
the base section is a cone. This is the shape utilized for numerous 
studies and is the shape of the slice pod. The generic figure is 


illustrated in Figure 3. 





Nose: Ellipise — Parallel Mid-body: Cylinder Base: cone shape 
°~L << |, <r 
Ln Lm Lb 
LB 
Lf 








Xm: Geometric Midpoint 
Leb: Center of Mass 





%: Individual section geometric center 
Figure 3. Geometric parameter description for body of revolution. 
The various methods utilized for the estimation of 


hydodynamic coefficients were compared by LT Wolkerstofer [2}. 
5 





The body of revolution hydrodynamic coefficients depend on 
semi-empirical relations to account for the viscous or vortex effects 
on the body. The methods based primarily on geometric 
considerations are applied in this parametric study. The 
hydrodynamic coefficients will be determined from the United States 
Air Force Data Compendium (USAF DATCOM) method which is 
summarized by Peterson [3]. The acceleration hydrodynamic 
coefficients will be determined by the methods of Humpreys and 
Watkinson [4]. 

To initiate the parametric study we needed to 
non-dimensionalize a number of characteristics of the body of 
revolution. The lengths of each section of the body are described by 


the fractional amount of the total length as in the nose fraction. 


| Ss (1) 


The slenderness ratio is defined using the total length. 


_ |e 
dd 


The volume the body of revolution would be determined from 


S (2) 


the sum of the three sections using standard volume formulas but 
with all the lengths referenced to the total length of the body via 


fractional amounts. 





2 
V= (a Jer. +3F,,+F,) (3) 








From this relation it is clear that the mid-body fraction is the 
dominant term in determination of the volume amount. This is 
consistent with the realization that the mid-body diameter is the 
maximum diameter for the body and constant mid-body length. The 
determination of the volume is important because it is the volume 
which determines the buoyant forces resulting from the shape. A 
larger volume will have a larger buoyant force. 

When the volume is non-dimensionalized the diameter of the 


body of revolution will be determined in the following relation. 


d=./(12V) /(nl,)(2F, +3F,, +F,) (4) 


The geometric center of the body of revolution will depend on 
the lengths of the individual sections. We utilized standard formulas 
for each section to find the geometric center for each section and 
then combined then to determine the offset from the geometric 


midpoint of the mid-body section. 


2 
a 21,(2+2*)—1,( 2+ | (5) 
12V 2° 8 2° 4 


This offset from the center of the mid-body section is then applied to 








find the geometric center for the body which is the same as the 


center of bouyancy for the body in horizontal motion. 


L,=l,+—-X, (6) 





The axial position where the flow becomes predominantly 
viscous is a function of the overall length determined from the point 
on the body of maximum slope [3]. The axial position, Ly is utilized 
for the determination of a few parameters affecting the 


hydrodynamic coefficients. 
1, =0.905 x1, (7) 


The geometry of the body at Ly is illustrated in Figure 4. 





Stb: Profile Area of Base 


Sv: Cross-section 
Area of Base 






2r 


LB -Lv 
Lv en ——_—_—————————e 


Figure 4. Sy and Stb illustration. 


The maximum profile area is determined at the mid-body 
section where the diameter is maximum. 
_ 1d* 


—— 8 
b A (3) 








The radius of the body at the axial position Ly is determined 


by realizing the base profile is a triangle. 


ies : 


The cross sectional area at the axial position is therefore 





S.=ar° (10) 


and the profile area of the base at the axial position is as follows: 
Sy =1(I5 -1,) (11) 


The drag coefficient for the body is initially established by 
Fidler and Smith [5] and utilized by Wolkerstofer [2]. The value of 


the drag coefficient will be modified within the parametric study in 


Chapter III. 
C, =0.29 (12) 
B. HYDRODYNAMIC COEFFICIENT CALCULATION 


From the basic physical description of the system the 
hydrodynamic parameters can now be determined. 

The Lamb's coefficients (kj and k2) of inertia are for a prolate 
ellipsoid in axial and cross flow. These parameters affect the 


lift/angle of attack curve slope [3]. 





k, —k 
C,, = 25,24 (13) 
b 


The hydrodynamic coefficient Yy is the normal force/angle of 


attack curve slope [3]. 
aS | 
opie aa lc, - C, (14) 
B 


The pitching moment/angle of attack curve slope was 
simplified using MATHCAD and its symbolic manipulator and is 


based on the cross sectional area of the body of revolution [3]. 





oa 2ke=h), 
Sle 
~—d? cee —I,) md? 
A a 
12 121? (15) 


3X yb, + 61,X,b, +215 — 311; — 31,1; 
$6X ql oly — OX qi -l, + 3U:1, —3x,,1¢ +1; 


The hydrodynamic coefficient Ny is the pitching moment 


coefficient [3]. 


N ime = Cay X= (16) 
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The lift/pitch rate curve slope is a function of the lift/angle of 


attack curve slope [3]. 


C,, =c, [1-24] (17) 


B 


The rotary derivative hydrodynamic coefficient Yy is the 


normal force/pitch rate coefficient [3]. 





C5 

) gemree i 4, (18) 

l 

The pitching moment/pitch rate curve slope [3]. 
2 
1-2] (“Hs a) 
l Sale 

C _ C B th’ B , (19) 





mq ~ —ma (1-%=]-[ V 
L Sale 


The rotary derivative hydrodynamic coefficient Ny is pitching 


moment/pitch rate coefficient [3]. 


a (20) 


The acceleration hydrodynamic coefficients Yy' and Ny’ are 


based on the work of Humphrys and Watkinson [4]. 
11 








Y = 


vdotprime ke ( 21 ) 
B 


N odotprime — = N detielins (22) 
The rotary acceleration hydrodynamic coefficient 1s modified 


by the mass moment of inertia of the displaced fluid about the z axis. 
lp 
2 
Ly = | pS(x)(x,, —x)° dx (23) 
0 


The rotary acceleration hydrodynamic coefficient N; is 


determined from the following relation [4]. 


2k, Lay 
pl; 


rdotprime ~ 


(24) 


C. PROGRAM VALIDATION 


We computed the values of the hydrodynamic coefficients 
using a MATLAB program and utilized the generalizations that will 
be required for the parametric study. The results were compared to 


the DATCOM SUBOFF data collected by LT Wolkerstofer [2] and the 


results are listed in Table 1. 








Function Error 


Hydrodynamic 
Coefficient 
Vv -Q. 
-0.0136 | -0.0136 
-0.0014 | -0.0014 


Po 
a mak 
0 
<< - 
ns 
Ls 









Table 1. DATCOM and program comparison. 





The results of the program agree with the data compiled from 


the actual model tested using the semi-empirical methods. The 


modifications to simplify the analysis based on the given shapes 


provided accurate results. The magnitude of the Ny’ error is due to 


the small value of the Nr' coefficient. 
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il. PARAMETRIC STUDY 


A. INTRODUCTION 


In conducting a parametric study of the body of revolution, we 
needed to decide which parameters to vary, which to hold constant 
and what effects these variations would have on the model. The 
primary goal was to find non-dimensional parameters that would 
define the shape sufficiently to accurately predict the hydrodynamic 
coefficients. The non-dimensional parameters would therefore not 
be constrained to a given particular size and shape. 

The data which completely describes the body of revolution 
include the sectional lengths and nominal (mid-body) diameter. 
From this the volume, cross sectional areas, and geometric centers 
could be calculated and the resulting hydrodynamic coefficients. 


This initial program was limited to specific non-variable inputs only 





which produced single case results. 

The initial steps in broadening the program were due to 
non-dimensionalizing the sectional lengths while still specifying the 
overall length of the body. The diameter was determined from the 
Slenderness ratio. In describing the body in this manner, we 
determined that the hydrodynamic coefficients were independent of 
the specific overall length of the body as the bodies parameters were 
all proportional to the length. We could therefore, alter the sectional 
fractional values to determine the effects on the hydrodynamic 


coefficients. 
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In order to study these variations in a controlled manner we 
decided to hold the volume of the body as a constant initially. With 
the volume as a constant, describing the sectional fractional values 
would determine the remaining parameters. In addition, only the 
nose and mid-body fractions would be altered since the base fraction 
would also be defined by default. 

From this basis we pursued two approaches, one involving a 
constant diameter and varying lengths and the second involving a 
constant length and varying diameter. This vectoral approach was 
later modified to allow varying lengths utilizing the non-dimensional 


volume/length? parameter. 


B. PHYSICAL REVIEW OF PARAMETER MODIFICATIONS 


In order to appreciate the multitude of possible parameter 
modifications Figures 5 through 8 are presented to solidify the 
physical changes incorporated in the parametric study. Figure 5 is 
the basic shape of the body of revolution. We have repeated this 


shape behind the modified shapes for Figures 6 through 8. 


Figure 5. Basic shape. 
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Figure 6. Lower Fp and slenderness ratio (shorter and fatter). 





Figure 8. Lower Fp, higher Fm, lower Fp, constant slenderness ratio. 


C. COMPUTER MODEL DEVELOPMENT 


In the process of the parametric study, the geometry of the 
body would be altered which would affect the assumptions and 
variables of the semi-empirical method utilized. The primary 
concern was with the drag coefficient. The changed shape will 


induce a different drag. We researched the significance and 
17 





magnitude of these modifications on the hydrodynamic coefficients. 
The value of drag coefficient was a constant 0.29 for various nose 
and base configurations in Fidler and Smith [5] but at the extreme 
sectional fractions it must be modified. We utilized the values from 
Herner [6] in the extreme cases. When the nose length to diameter 
ratio was less than 1 the coefficient of drag was smoothly increased 
to the maximum value of 0.82. When the base length to diameter 
ratio was less than 1 the coefficient of drag was smoothly increased 
to the maximum value of 0.64. When incorporated into the program 
this parameter variation had minimal effect on the final result since 
the other components prescribed by the parameter values dominated 
the solution trend and our specific area of interest 1s not at the 
sectional fractional extreme. 

The data format illustrated in Figure 9 demonstrates the limits 
of the model and is useful in interpreting the mesh graphs. 

Since the nose, mid-body and base fractions together comprise 
100% of the body the values in the lower right hand side diagonal 
are not possible combinations and an arbitrary constant value is 
utilized to complete the matrix and provide a visible floor to the 
mesh graphs. Appendix A is the MATLAB program that generated 


the matrix data with the full range of nose and mid-body fractons. 
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Figure 9. Matrix interpretation for data analysis. 
D. HYDRODYNAMIC COEFFICIENT MESH GRAPHS 


The Figures 10 through 17 were generated with an automatic 
scaling feature. This allows for a rough comparison of the magnitude 
of each parameter to be compared to another parameter. The Yy' 
variations are significant while the acceleration hydrodynamic 
coefficients variations in magnitude are relatively small. The graphs 
of Ny' and Ny' included a clipping feature to allow for the three 
dimensional viewing. The clipping only involves truncating the 
values after the solution begins to reach an extreme value. 

The Figures 19 through 32 are scaled to determine the shape of 
the hydrodynamic coefficient and its variations over the entire range 


of sectional fractions. Each hydrodynamic coefficient mesh graph is 
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coupled with a two dimensional graph which is not clipped and 
illustrates the varying parameter in a simpler and therefore more 


clean environment. 
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Figure 11. Yvprime standard scale mesh graph. 
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Figure 12. Nvprime standard scale mesh graph. 
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Figure 14. Nrprime standard scale mesh graph. 
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Figure 16. Nvdprime standard scale mesh graph. 
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Figure 19. 


Zz9 





PD 


i 2 


a i 2 


a 


i FE. he a 
8 Sa Nek 
: iS 


a oF 


Pe i 


pee ew we we eee th ww tk tt Ly wae we ed 


) Tp) st 
=“ +P) = 
© ™~ © 
ro) ©o von) 
| oO 1 
SWUdAN 








Mid-body Fraction 


0 


Nose Fraction 


up scale mesh graph. 


Nvprime close- 
-0.013 


Figure 20. 


Nose Fraction (%) 





Te) 
o 
se 
° 
? 


OWUdAN 





-0.0145 


0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Mid-body Fraction 


0.1 


Nvprime variations with mid-body fraction. 
30 


Figure 21. 
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Figure 23. Yrprime variations with mid-body fraction. 
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Figure 24. 
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Nrprime variations with mid-body fraction. 
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Figure 25. 
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Figure 26. Yvdprime close-up scale mesh graph. 
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Figure 27. Yvdprime variations with mid-body fraction. 
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Figure 28. Nvdprime close-up scale mesh graph. 
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Figure 29. Nvdprime variations with mid-body fraction. 


34 








x 10 


Nrdotprime 





Nose Fraction Mid-body Fraction 


Figure 30. Nrdprime close-up scale mesh graph. 
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Figure 31. Nrdprime variations with mid-body fraction. 
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E. VOLUME AND LENGTH VARIATIONS 


The mesh graphs of the hydrodynamic coefficients were 
generated holding the volume of the body of revolution constant. We 
now explored variations in the volume and length parameters of the 
body of revolution to see the effects on the hydrodynamic 


coefficients as illustrated in Figure 32. 


Figure 32. Variations of V/L? ratio drawing. 


We took the basic body and coupled the volume and length to 
determine if a non-dimensional volume/length? ratio was valid. The 
non-dimensional results matched the original results. The value of 
volume/length? for the body was 8.023x10-3. With this value as an 
average value we varied the volume/length> ratio from 
6 to 10x10-3. The range is noted on all figures as 6 to 10 for 
simplicity. The smaller value is due to a smaller volume or a larger 
overall length. Figures 33 through 46 investigate these variations at 
nose fractions of 10% and 20%. For the vast majority of the entire 
range variations of volume/length> result in linear trends in all 


CaSes. 
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Figure 33. Yvprime V/L3 variations at 10% nose fraction. 
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Figure 34. Yvprime V/L° variations at 20% nose fraction. 
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Figure 35. Nvprime V/L2 variations at 10% nose fraction. 
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Figure 36. Nvprime V/L3 variations at 20% nose fraction. 
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Figure 37. Yrprime V/L> variations at 10% nose fraction. 
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Figure 38. Yrprime V/L? variations at 20% nose fraction. 
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Figure 39. Nrprime V/L> variations at 10% nose fraction. 
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Figure 40. Nrprime V/L3 variations at 20% nose fraction. 
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Figure 41. Yvdprime V/L> variations at 10% nose fraction. 
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Figure 42. Yvdprime V/L? variations at 20% nose fraction. 
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Figure 43. Nvdprime V/L3 variations at 10% nose fraction. 
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Figure 44. Nvdprime V/L> variations at 20% nose fraction. 
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Figure 45. Nrdprime V/L> variations at 10% nose fraction. 
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Figure 46. Nrdprime V/L3 variations at 20% nose fraction. 
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IV. DETERMINATION OF FUNCTIONAL COEFFICIENTS 


A. INTRODUCTION 


The hydrodynamic coefficients for various combinations of 
sectional fractions and volume/length ratios have now been 
computed and graphed. The surface functions were generated using 


the semi-empirical methods under the following general formula. 


He = FF, -Fa | | (25) 


The goal in predicting the hydrodynamic coefficients would be 
achieved if we could determine the relationship along the surface 
function based on the variable parameters. 

In order to minimize error and simplify the functional 
relationships we limited the range of the variable parameters and in 
each case the range of the parameters covers the majority of the 
cases of interest. The nose fraction range is 5% to 25%. The 
mid-body fraction is 40% to 60%. The volume/length? range is 
6 to 10x10-3. 

In reviewing Figures 18-31, in our range of interest, the curves 
appear to correspond to a 2nd degree polynomial. In reviewing 
Figures 33-46, in our range of interest, the curves have a definite 
linear relationship. The difficulty would be in determining the 


surface functional relationship based on two parameters since 


45 





currently there is no built-in subroutine in the MATLAB program 


toolbox in use to solve this problem. 
B. PREDICTION EQUATION 


Hard copy updates and distribution of improvements in 
software always involves time delays. In corresponding with the 
Mathworks company I learned that a collection of m-files existed on 
the internet. I utilized some of these m-files in the development of 
my graphs (the legend command is not included on the Macintosh 
MATLAB professional version). I downloaded a surface equation 
least squares curve fitting program that I could modify to determine 
the coefficients of the surface function. I tested the program out 
using a number of different surface functions and concluded the 
program accurately predicted the correct polynomial function. 
Appendix B includes a sample program. 

With the functional coefficients determined the hydrodynamic 


coefficient prediction equation would be as follows: 
HC=A,F°+A,F.F.+4A,F° +A,F,+A,F,, +A, 


V (26) 
ral a-G) 


The coefficients for the equations are listed in Table 2 for the 


nose, mid-body fraction and volume/length? ratio. The constant C1 


is 8.023x10-3 and is the nominal value for volume/length? ratio. 
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Figures 47 through 60 compare the theoretical surface 
functions to the predicted surface functions using the same axial 
scaling. Figures 60 through 67 are included to demonstrate the 
percentage error of the prediction equation from the theoretical 
surface function. In general, the percentage error is small with the 
exception of the Yr percentage error but that is due to the relatively 
smaller magnitude of the Yr hydrodynamic coefficient. The 


acceleration hydrodynamic prediction equations are very accurate. 





0.0002 | 0.0007 | 0.0007 | -0.0008 | -0.0016 |} -0.0144 | -1.8067 
-0.0002 | -0.0007 | -0.0007 | 0.0008 | 0.0016 | 0.0144 | 1.8067 
-0.0031 | -0.0046 | -0.0021 | 0.0031 | 0.0024 | -0.0013 | -0.0808 





Table 2. Functional Coefficients for prediction equation. 


The coefficients of Table 2 are calculated for the following 
nominal values: nose fraction 15%, mid-body fraction 50%, and the 


volume/length3 ratio 8.023x10-3. These parameters represent the 
mid-range values. 
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Figure 47. Yvprime theoretical mesh graph. 
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Figure 48. Yvprime predicted mesh graph. 
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Figure 50. 


49 














oe 
ot * 





x10~ 
9) 
es 
1 
: ORE SS Os 
TS 
» 
E 
a. 
= 3 : 
Baal ceeere : 
5 Oe oe 
0.05 
0.1 
0.15 
ws bates an 0.55 
0.2 2 0.5 
0.45 
0.25 4 
Nose Fraction Mid-body Fraction 


Figure 51. Yrprime theoretical mesh graph. 
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Figure 52. Yrprime predicted mesh grpah. 
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Figure 54. Nrprime predicted mesh graph. 
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Figure 57. 


Mid-body Fraction 


53 


Nose Fraction 


Nvdprime predicted mesh graph. 


Figure 58. 





Nrdprime 


Figure 59. 


Nrdprime Calculation 


Figure 60. 


Nrdprime predicted 


mesh graph. 


54 


~4 ‘aot 
x 10 at 
; ae EES 
Se 
FEF! 
NE 
LEELA DPA 
SS SSRN EMAL 
ae SN UN <S a 
- SSSR 
SSS SSS 
Laer AY A TARY 
eo . “SSS 
. Os ae “S SVAN ‘ . 
tne i oo Ca \ a 
Oe en Se eT = OA RA 
: Sy cAh ok 
a I EP ; Sk 
O41. eng ee. 
0.15 ae 
ame ee 0.55 
0.2 0.5 
0.45 
0.25 4 . 
Nose Fraction Mid-body Fraction 
Nrdprime theoretical mesh graph. 
nn ie 
x 10 cae : 
i ees SSESE_R 
SSS << 
LEE 
EES SSSA 
LLP 
LAZZARO, 
SSSR SSS YAS ee 
Se SOS LSS Ae 
eee RESELLERS AN 
cakes Seas SSN SS a 
ey : el as oe : A SC 
; . ae ae ~ a 
Bc ase BT ete ™s a \ aN 
pet a cae a ASA, 
¥ “. : es: 
OS ee OR Ue 
0.15 
0.55 
0.2 0.5 
0.45 
0.25 94 
Nose Fraction Mid-body Fraction 


0.6 


0.6 





0.55 





c 
2 
ro) 
a] 
S ‘ 
: =" 4 
in Ss : : 
| pa 
ro) gs : : : 
+ pa : 
RS @ 
= = ; 
; OD 
LO ‘is : 
v SS : 
: : oD) Ww) 
~~ ; : 
3 o 
= a, 
3 ro) 
; S 
So o. ; : 
S a g 
Lf) oO : 
N iS 
o>) = 
; re) 
iS a § . 
oO - @ 
a @ s E 
: @ 
Te) © = : = 
banal Tm — : : 
° a) hy E 
8 jar 
3 
Ye) 
oO oO © 2 
g eo 2 38 


JOIIF JUGS SWHCGAN 


Figure 61. 
Figure 62. 
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Figure 63. Yrprime equation percentage error. 
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Figure 64. Nrprime equation percentage error. 
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Figure 65. 


Nvdotprime equation percentage error. 


Figure 66. 
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Figure 67. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 
1. Summary Review of Parameter Variations 


The matrix variations of the nose and mid-body fractions 
resulted in a mixture of surface functions for the hydrodynamic 
coefficients. There was not an overriding trend that was consistent 
with all the hydrodynamic coefficients. The trends were all 2nd 
order polynomials for sectional fractions and linear for 
volume/length ratio variations. The larger prediction error for the 
Yr is due to the relatively smaller magnitude variation in the 
hydrodynamic coefficient. The acceleration hydrodynamic 
coefficients were well behaved functions and the prediction 
equations are highly accurate. The acceleration hydrodynamic 
coefficients also did not vary greatly in magnitude. Therefore, the 
designer can alter the sectional fractions without too great a concern 
for the effects due to acceleration. 

A lower volume/length value reduces the magnitude of the 
hydrodynamic coefficient and this is consistent for all the cases 
studied. For a constant volume, a lower volume/length implies a 
longer slimmer body, and a higher slenderness ratio. This would 
imply the effect of the diameter on the body has a greater impact 


than the length. 


59 








2. Final Comparison 


We compared the program results to the expected parameters 
of the SUBOFF Body and the Slice Pod The results are listed in 
Table 3. 


Slice Pod 
¥y | -0-0058 | -0.0054 | 6.8 | -0.0141 | -0.0153 | 7.8 


Table 3. Program Comparison. 














Overall, the comparison is acceptable with the understanding 
that the percentage error with Ny is due to its relatively small value. 


The Yr percentage error is within the tolerance of the program. 
B. RECOMMENDATIONS 
The parametric study of the body of revolution while singular 


in scope presents the ability to pursue variations in the design of 


components and their effect on the hydrodynamic coefficients. 
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Recommendations for further research in this area are as follows: 


¢ Modify the program to evaluate different shapes. 
¢ Evaluate motion in the vertical plane. 
¢ Calculate stability criteria and general maneuvering 


performance for the surface mesh functions. 
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APPENDIX A. MATLAB ROUTINE FOR DETERMINING 


HYDRODYNAMIC COEFFICINETS WITH VARYING NOSE 


A. 


AND MID-BODY FRACTIONS 


VDCURVES.M PROGRAM 


% LCDR Eric Holmes 


% 
% 


Constant: Volume/Length’3 ratio 
Thesis 


clear 
clear global 


%o 


Basic Parameters 


IB = Total body length 

In = Nose length 

Im = Middle body length 

lb = Base body length 

d = diameter 

If = forebody length (lb + lm) 

Fn = Nose fractional length of total body length 

Fm = Mid-body fractional length of total body length 

r_Fm = maximum fraction of mid-body 

Fb = Base fractional length of total body length 

N,Lk,t = index counters 

Sb = maximum cross sectional area of body (assumimg diameter 

d and length IB) 

V = Volume of the body 

xm = Geometric middle of body 

xcb = Geometric offset from the center of the middle body 
for gravity 


Icb = Geometric center of gravity 
rho = density of water 
e = Munk coefficient 


Bo = Munk coefficient 
Ao = Munk coefficient 


k1 = Lamb's inertial coefficient 
k2 = Lamb's inertial coefficient 
kb = Lamb's inertial coefficient 


Cdo = Drag coefficient at zero angle of attack 
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% lv = length where viscous flow dominates 

% xr = radious at length lv 

% R = range for graphs 

% S = Slenderness ratio matrix 

% Sratio = slenderness ratio 

% Sv = cross sectional area at length lv 

% Stb = platform sectional area at lenght lv 

% Cla = Lift/Angle of attack curve slope 

% Cma = Pitching moment/Angle of attack curve slope 
% Cmq = Pitching moment/pitch rate curve slope 

% Clq = Lift/Pitch rate curve slope 

% Iydf = mass moment of inertia of displaced fluid mass 
% Tydf_In = nose body component of Iydf 

% lydf_lIm = middle body component of Iydf 

% lydf_lb = base body component of Iydf 


% Hydrodynamic coefficients 

% Yvprime = normal force coefficient 

% Nvprime = pitching moment coefficient 

% Yrprime = normal force/pitch rate coefficient 

% Nrprime = pitching moment/pitch rate coefficient 

% Zwdotprime = acceleration coefficient (axisymetric to Yvdotprime) 
% Yvdotprime = acceleration coefficient along y axis 

% Nvdotprime = acceleration coefficient causing yawing moments 

% Nrdotprime = acceleration coefficient in addition to Iz 


% Yvp,Nvp,Yrp,Nrp,Yvdp,Nvdp,Nrdp are hydrodynamic 

% coefficient matrices each combining the values of their respective 
% coefficient: Yvpprime => Yvp 

% Yvp_d,Nvp_d,Yrp_d,Nrp_d,Yvdp_dNvp_d,Nrp_d are hydrodynamic 
% coefficient matrices maintained as raw data. 

% (The matrix filler is -1.0) 

% Yvp_s,Nvp_s,Yrp_s,Nrp_s,Yvdp_s,Nvp_s,Nrp_s are hydrodynamic 

% coefficient matrices that have been clipped for graphical 

% presentation 


% Initializing data and empty matrices 
S =[]; Yvp = []; Nvp = (]; Yrp = (J; Nrp = []; Yvdp =[]; Nvdp = []; Nrdp = []; 
R = (0:0.01:1); 
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Ve 
IB = 4.9941; 
rho = 62.4/32.174; 


% Drag Function Development 


In_d = [0 .5 1 1.5 2 2.5 3]; 
cd_Ind = [0.82 .45 .29 .29 .29 .29 .29]; 
cln = polyfit(in_d,cd_Ind,4); 


Ib_d = [0 511.52 2.5 3]; 
cd_Ibd = [0.64 .38 .29 .29 .29 .29 .29}; 
clb = polyfit(Ib_d,cd_lbd,4); 


% Establishing fractional values for the body 
Fn = (0:0.01:1); 

for k = 1:101; 

r_Fm = 1 - Fn(k); 


Fm = (0:0.01:r_Fm); 
Fb = 1 - Fn(k) - Fm; 


% Diameter calculation 

d = ((12*V)./((pi*1B).*(2*Fn(k) + 3*Fm + Fb))).40.5; 
% Length calculation 

In = Fn(k).*IB; 

lm = Fm.* IB; 

lb = Fb.*1B; 

If = In + Im; 


% Slenderness ratio calculation 


Sratio = 1B./d; 
S = [S Sratio zeros(1,k-1)]; 


% Body calculation 
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Sb = (pi.*d.42)/4; 
xm = [B/2; 


global d IB If In Im 1b xm rho I 

% Center of gravity and bouyancy 

xcb = (pi*d.42./12).*(2.*In.*(m/2 + 3*In/8) - 1b.*Um/2 + 3*1b/4))/V; 
Icb = In + Im/2 - xcb; 

% Munk coefficients 


e = (2./IB).*(((IB.42)/4 - Sb/pi)).40.5; 

Bo = (1./e.42) - ((1 - €.42)./(2*e.43)).*log((1+e)./(1-e)); 

Ao = ((1-(e.42))./(e.43)).*(log((1+e)./(1-e)) - 2*e); 

k1 = Ao./(2 - Ao); 

k2 = Bo./(2 - Bo); 

kb = ((e.44).*(Bo - Ao))./((2 - €.42).*(2*%e.42 - (2 - e.42).*(Bo - Ao))); 


% Normal force coefficients 
% Drag coefficient determination 


for 1 = 1:102 - k 
if In/d@) < 1 
Cdo(i) = cln(1)*(In/dq))*4 + cln(2)*(n/d@i))43 + cln(3)*(in/d(i))42 + ... 
cln(4)*(In/d@)) + cln(5); 
else 
Cdoi) = 0.29; 
end 
if lbG@)/di) < 1 
Cdod) = clb(1)*(b(i)/d(@i))*4 + clb(2)*(db(i)/dqG))43 + =... 
clb(3)*(Ib(i)/dq))*2 + clb(4)*(bG)/d(@i)) + clb(5); 
else 
Cdo(1) = 0.29; 
end 
end 


lv = 0.905*IB; 


r = (d/2).*(1 - (lv - If)./1b): 
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SV = pi.*r.42; 
Cla = (2*(k2 - k1).*Sv)./Sb; 
Yvprime = (-Sb./IB.42).*(Cla + Cdo); 


Yvp = [Yvp Yvprime -(ones(1,k-1))]; 

Cma = (2*(k2-k1)./(Sb*1B)).*((-pi.*d.42./12).*(3*xm - In) + ... 
((p1.*d.42)./(12*1b.42)).* ... 
(-3*xm.*lv.A2 + 6*lv.*xm.*lb + 2*lv.43 - 3*1f.*lv.42 - 3*lb.*lv.A2 + ... 
6*lv.*xm.*lf - 6*1f.*xm.*Ib + 3*]b.*1f.42 - 3*xm.*1f.42 + If.43)); 

Nvprime = (Sb.*Cma)./IB.42; 

Nvp = [Nvp Nvprime -(ones(1,(k-1)))]; 

% Rotary force coefficients 


Stb = r.*(1B - Iv); 


Cmgq = Cma.*((1 - xm./IB).42 - V*(Icb - xm)./(Stb.*1B.%2)),/... 
(1 - xm./IB) - (V./(Stb.*1B))); 


Nrprime = -(Sb.*Cmq)./IB.‘2; 


Nrp = [Nrp Nrprime -(ones(1,(k-1)))]; 
clear Nrprime 


Clq = Cla.*(1- xm./IB); 

Yrprime = (-Sb.*Clq)./IB.%2; 

Yrp = [Yrp Yrprime -(ones(1i,(k-1)))]; 
% Acceleration coefficients 
Zwdotprime = (2*k2*V)./(1B.43); 
Yvdotprime = -Zwdotprime; 


Yvdp = [Yvdp Yvdotprime -(ones(1,(k-1)))]; 
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for I = 1:102-k; 
Iydf_Indl) = quad8('vdn',0,1n); 
lydf_Im(D= quad8(‘vdm',In,1f(1)); 
Iydf_Ib(1) = quad8('vdb',lf(D,1B); 
end 
Iydf = Iydf_In + lydf_lm + Iydf_lb; 
Nrdotprime = (-2.*kb.*Iydf)./(rho.*1B.45); 
Nrdp = [Nrdp Nrdotprime -(ones(1,(k-1)))]; 
N = isnan(Nrdp); 
I = find(N > 0); 
for t = 1:length(I) 
Nrdp(i(t)) = -1.0; 


end 


clear d lf In lm lb xm I 

clear Iydf_In Iydf_lm lIydf_lb lIydf 
clear global 

IB = 4.9941; 

tho = 62.4/32.174: 


end 


Nvdp = -Yvdp; 
Nvdp = reshape(Nvdp,101,101); 


= reshape(Yvp,101,101); 
Nvp_d = reshape(Nvp,101,101); 


Yrp_d reshape(Yrp,101,101); 
Nrp_d reshape(Nrp, 101,101); 
Yvdp_d reshape(Y vdp, 101,101); 


Nvdp_d = reshape(Nvdp,101,101); 
= reshape(Nrdp,101,101); 


save Yvp_d Nvp_d Yrp_d Nrp_d Yvdp_d Nvdp_d Nrdp_d 
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% Format for presentation (all analysis complete) 
Nvp = reshape(Nvp,101,101); 


figure(1) 

plot(R,Nvp(1:101,6),'c', R,Nvp(1:101,11),'r+', R,Nvp(1:101,16),'b*', ... 
R,Nvp(1:101,21),'go') 

xlabel('Mid-body Fraction’) 

ylabel(’Nvprime') 

text(0.1,-0.0132,'Nose Fraction (%)') 

legend(' 5','10','15','20') 

axis({O 1 -0.0145 -0.013]) 

hold on 

plot(R,Nvp(1:101,11),'r:', R,Nvp(1:101,16),'b:'", R,Nvp(1:101,21),'g:') 

hold off 

pause 

print -depsc2 nvp 


Nvp = reshape(Nvp,1,10201); 

k = find(Nvp > -0.01); 

Nvp(k) = -0.02*(ones(length(k),1)); 
Nvp = reshape(Nvp,101,101); 


k = find(Nrp < -0.01); 
Nrp(k) = ones(length(k),1); 
Nrp = reshape(Nrp,101,101); 


figure(2) : 

plot(R,Nrp(1:101,6),'c’, R,Nrp(1:101,11),'r+', R,Nrp(1:101,16),'b*', ... 
R,Nrp(i:101,21),'go') 

xlabel('Mid-body Fraction‘) 

ylabel('Nrprime') | 

text(0.2,0.0025,'Nose Fraction (%)') 

legend(' 5','10','15','20') 

axis({O 1 -0.0025 0.003]) 

hold on 

plot(R,Nrp(1:101,11),'r:'", R,Nrp(1:101,16),’b:', R,Nrp(1:101,21),'g:') 

hold off 

pause 

print -depsc2 nrp 
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Nrp = reshape(Nrp,1,10201); 

k = find(Nrp > 0.05); 

Nrp(k) = -0.02*(ones(length(k),1)); 
Nrp = reshape(Nrp,101,101); 


figure(3) 

plot(R,Nrdp(1:101,6),'c’, R,Nrdp(1:101,11),'r+', R,Nrdp(1:101,16),'b*', ... 
R,Nrdp(1:101,21),'go') 

xlabel(Mid-body Fraction’) 

ylabel(’'Nrdotprime') 

text(0.6,-0.00015,'Nose Fraction (%)') 

legend(’ 5','10','15','20') 

axis([O 1 -0.0012 0]) 

hold on 

plot(R,Nrdp(1:101,11),'r:',, R,Nrdp(1:101,16),'b:', R,Nrdp(1:101,21),'g:’) 

hold off 

pause 

print -depsc2 nrdp 


k = find(Nrdp < -0.02); 
Nrdp(k) = -0.02*(ones(length(k),1)); 
Nrdp = reshape(Nrdp,101,101); 


k = find(Yvp<-0.02); 
Yvp(k) = -0.02*(ones(length(k),1)); 
Yvp = reshape(Yvp,101,101); 


k = find(Yrp < -0.02); 
Yrp(k) = -0.02*(ones(length(k),1)); 
Yrp = reshape(Yrp,101,101); 


k = find(Yvdp < -0.02); 
Yvdp(k) = -0.02*(ones(length(k),1)); 
Yvdp = reshape(Yvdp,101,101); 


figure(4) 
plot(R, Yvp(1:101,6),'c', R,Yvp(1:101,11),'r+', R,Yvp(1:101,16),'b*', ... 
R, Yvp(1:101,21),'go') 
xlabel(/Mid-body Fraction’) 
ylabel( Y vprime') 
text(0.15,-0.0075,'Nose Fraction (%)') 
legend(' 5','10','15','20') 
70 





axis({0 1 -0.01 -0.004]) 

hold on 

plot(R, Yvp(1:101,11),'r:', R,Yvp(1:101,16),'b:', R,Yvp(1:101,21),'g:') 
hold off 

pause 

print -depsc2 yvp 


figure(5) 

plot(R, Yrp(1:101,6),'c', R,Yrp(1:101,11),'r+', R, Yrp(1:101,16),’b*', ... 
R, Yrp(1:101,21),'go') 

xlabel('Mid-body Fraction‘) 

ylabel('Yrprime') 

text(0.15,-0.0025,'Nose Fraction (%)') 

legend(' 5','10','15','20') 

axis({0 1 -0.005 0]) 

hold on 

plot(R, Yrp(1:101,11),'r:', R,Yrp(1:101,16),'b:'", R,Yrp(1:101,21),'g:’) 

hold off 

pause 

print -depsc2 yrp 


figure(6) | 

plot(R, Yvdp(1:101,6),’c', R,Yvdp(1:101,11),'r+', R,Yvdp(1:101,16),'b*', ... 
R, Yvdp(1:101,21),'go') 

xlabel(‘Mid-body Fraction’) 

ylabel('Y vdotprime') 

text(0.5,-0.0146,'Nose Fraction (%)') 

legend(' 5','10','15','20') 

axis({0 1 -0.0154 -0.0145]) 

hold on 

plot(R, Yvdp(1:101,11),'r:', R,Yvdp(1:101,16),'b:', R,Yvdp(1:101,21),'g:') 

hold off 

pause 

print -depsc2 yvdp 


figure(7) : 
plot(R,Nvdp(1:101,6),'c’, R,Nvdp(1:101,11),'r+', R,Nvdp(1:101,16),... 
'b*', R,Nvdp(1:101,21),'go') 
xlabel('Mid-body Fraction’) 
ylabel('Nvdotprime') 
text(0.4,0.015,'Nose Fraction (%)') 
legend(’ 5','10','15','20') 
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axis([O 1 0.0145 0.0154]) 

hold on 

plot(R,Nvdp(1:101,11),'r:', R,Nvdp(1:101,16),'b:', R,Nvdp(1:101,21),'g:’) 
hold off 

pause 

print -depsc2 nvdp 


% Range change for mesh graphics 
R = (0:0.02:1); 


figure(8) 

S = reshape(S,101,101); 
mesh(R,R,S(1:2:101,1:2:101)), grid 
xlabel(‘Nose Fraction’) 
ylabel('Mid-body Fraction’) 
Zlabel('Slenderness Ratio’) 
view(60,30) 

print -depsc2 slender 


figure(9) 
mesh(R,R,Yvp(1:2:101,1:2:101)), grid 
Xlabel('Nose Fraction’) 
ylabel('Mid-body Fraction’) 
Zlabel('Y vprime’) 

view(60,30) 

print -depsc2 yvpm 


figure(10) 
mesh(R,R,Nvp(1:2:101,1:2:101)), grid 
xlabel('Nose Fraction’) 
ylabel('Mid-body Fraction’) 
Zlabel('Nvprime') | 
view(40,40) 

print -depsc2 nvpm 


figure(11) 
mesh(R,R, Yrp(1:2:101,1:2:101)), grid 
xlabel('Nose Fraction’) 
ylabel('Mid-body Fraction’) 
Zlabel('Yrprime') 
view(60,30) 
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print -depsc2 yrpm 


figure(12) 
mesh(R,R,Nrp(1:2:101,1:2:101)), grid 
Xlabel('Nose Fraction’) 
ylabel(Mid-body Fraction’) 
Zlabel('Nrprime') 

view(40,30) 

print -depsc2 nrpm 


figure(13) 

mesh(R,R, Yvdp(1:2:101,1:2:101)), grid 
Xlabel(‘Nose Fraction’) 
ylabel(’‘Mid-body Fraction’) 
Zlabel(' Y vdotprime’) 

view(60,30) 

print -depsc2 yvdpm 


figure(14) 
mesh(R,R,Nvdp(1:2:101,1:2:101)), grid 
Xlabel('Nose Fraction’) 
ylabel('Mid-body Fraction’) 
Zlabel(’Nvdotprime’) 

view(20,20) 

print -depsc2 nvdpm 


figure(15) 
mesh(R,R,Nrdp(1:2:101,1:2:101)), grid 
Xlabel('Nose Fraction’) 
ylabel(’Mid-body Fraction’) 

Zlabel( Nrdotprime') 

view(60,50) 

print -depsc2 nrdpm 


% Clipping of data for graphical presentation 


Yvp_s = reshape(Yvp,1,10201); 

k = find(Yvp_s<-0.01); 

Yvp_s(k) = -0.01*(ones(length(k),1)); 
Yvp_s = reshape(Yvp_s,101,101); 


Nvp_s = reshape(Nvp,1,10201); 
fae: 





k = find(Nvp_s > -0.013); 

Nvp_s(k) = -0.0145*(ones(length(k),1)); 
k = find(Nvp_s < -0.0145); 

Nvp_s(k) = -0.0145*(ones(length(k),1)); 
Nvp_s = reshape(Nvp_s,101,101); 


Yrp_s = reshape(Yrp,1,10201); 

k = find(Yrp_s < -0.005); 

Yrp_s(k) = -0.005*(ones(length(k),1)); 
Yrp_s = reshape(Yrp_s,101,101); 


Nrp_s = reshape(Nrp,1,10201); 

k = find(Nrp_s < -0.005); 

Nrp_s(k) = -0.005*(ones(length(k),1)); 
Nrp_s = reshape(Nrp_s,101,101); 


Yvdp_s = reshape(Yvdp,1,10201); 

k = find(Yvdp_s < -0.0154); _ 

Yvdp_s(k) = -0.0154*(ones(length(k),1)); 
Yvdp_s = reshape(Yvdp_s,101,101); 


Nvdp_s = -Yvdp_s; 


Nrdp_s = reshape(Nrdp,1,10201); 

k = find(Nrdp_s < -0.00125); 

Nrdp_s(k) = -0.00125*(ones(length(k),1)); 
Nrdp_s = reshape(Nrdp_s,101,101); 


figure(16) 

mesh(R,R, ¥vp_s(1:2:101,1:2:101)), grid 
Xlabel('Nose Fraction’) 
ylabel('Mid-body Fraction’) 

Zlabel('Y vprime’) 

view(60,30) 

print -depsc2 yvps 


figure(17) 
mesh(R,R,Nvp_s(1:2:101,1:2:101)), grid 
xlabel(‘Nose Fraction’) 
ylabel( Mid-body Fraction’) 
Zlabel(’Nvprime') 
view(40,40) 
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print -depsc2 nvps 


figure(18) 

mesh(R,R, Yrp_s(1:2:101,1:2:101)), grid 
Xlabel(‘Nose Fraction’) 
ylabel(’Mid-body Fraction’) 
Zlabel('Yrprime') 

view(60,30) 

print -depsc2 yrps 


figure(19) 
mesh(R,R,Nrp_s(1:2:101,1:2:101)), grid 
Xlabel('Nose Fraction’) 
ylabel('Mid-body Fraction’) 
Zlabel('Nrprime') 

view(40,40) 

print -depsc2 nrps 


figure(20) 
mesh(R,R,Yvdp_s(1:2:101,1:2:101)), grid 
Xxlabel('Nose Fraction’) 
ylabel(‘Mid-body Fraction’) 

Zlabel('Y vdotprime' ) 

view(60,30) 

print -depsc2 yvdps 


figure(21) 
mesh(R,R,Nvdp_s(1:2:101,1:2:101)), grid 
xlabel(‘Nose Fraction’) 

ylabel('Mid-body Fraction’) 
Zlabel('Nvdotprime’) 

view(40,30) 

print -depsc2 nvdps 


figure(22) 
mesh(R,R,Nrdp_s(1:2:101,1:2:101)), grid 
Xlabel(‘Nose Fraction’) 
ylabel(‘Mid-body Fraction’) 
Zlabel('Nrdotprime’) 
axis({0 1 0 1 -0.00125 -0.0005]) 
view(60,50) 
print -depsc2 nrdps 
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B. VDN.M PROGRAM 


% LCDR Eric Holmes 
% Thesis 
Ylydf_In Functional Integration 


% d = diameter 

% 1B = total body length 

% \n = nose length 

% \m = middle length 

% Ib = base length 

% Xm = geometric middle of the body 
% rho = density of water 


function I1 = vdn(x) 
global d IB If In lm lb xm rho I 


I1 = (Crho*pi.*d(1).42.*x.*(xm - x).42./(4.*1n‘2)).*(x - 2.*In)); 


C. VDM.M PROGRAM 


%LCDR Eric Holmes 
% Thesis 
%lydfi_Im Functional Integration 


% d = diameter 

% 1B = total body length 

% \n = nose length 

% \m = middle length 

% \b = base length 

% xm = geometric middle of the body 
% rho = density of water 


function [12 = vdm(x) 
global d IB If In lm Ib xm rho I 
I2 = (rho*pi.*d(1).42/4).*(xm - x).42; 
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D. VDB.M PROGRAM 


%LCDR Eric Holmes 
% Thesis 
Ylydf_lb Functional Integration 


% d = diameter 

% 1B = total body length 

% \n = nose length 

% Im = middle length 

% \b = base length 

% xm = geometric middle of the body 
% tho = density of water 

% \f = fore body (ib + Im) 


function [3 = vdb(x) 
global d 1B lf In Im lb xm rho I 


I3 = (pi*rho.*d(1).42.*(xm - x).42).*(b(I) - x + 1f(D)).42./(4*Ib(1).42); 


VW 


78 





APPENDIX B. SAMPLE MATLAB ROUTINE FOR DETERMINING 
FUNCTIONAL COEFFICIENTS FROM A SURFACE PROFILE 


A. SURF.M PROGRAM 


% LCDR Eric P. Holmes 
% Surface Equation Solver 


i= 2: 

x=[) y=0); 

load hcdata 

Yvp = Yvp_d(6:26,41:61); 
Yvp = Yvp(:); 


% Range 


xr = (5:25); 

for 1 = 1:21 
X = [x;xr]; 

end 

x = 0.01*x; 

X = x(*); 


yr = (40:60); 

for 1 = 1:21 
y = ly yr; 

end 

y = 0.01*y; 

y = y(); 


% Program operation 


n=n+1; 
ke 
A = zeros(size(x)); 
for i = n:-1:1, 
forj= 11 
ACK) = ((x.A0-j)).*(y.4G-1))); 
k = k+l; 
end 
end 
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p = (A\Yvp).' 

% Error calculation 

xX = reshape(x,21,21); 

y = reshape(y,21,21); 
Yvp = reshape(Yvp,21,21); 


for 1 = 1:21 
for } = 1:21 





Yvpeal(i,j) = p(1)*x(@.j)*2 + p(2)*x(ij)*yG@p + p(3)*yG,j*2 + .. 


p(4)*x,j) + p(S)*yGj) + p(6); 
end 
end 


for t=: 1:2) 
for j = 1:21 

perr(i,j) = (Yvpeal(i,j) - Yvpi,j))/YvpGi.j): 
end 

end 


perr = 100*perr; 


xp = x(1,:); 
yp = y(,1); 
figure(1) 


mesh(xp,yp,Yvp), grid 
Xlabel(Nose Fraction’) 
ylabel('Mid-body Fraction’) 
Zlabel('Yvprime') 
view(60,30) 

print -depsc2 yvp 


figure(2) 
mesh(xp,yp,perr), grid 
xlabel(Nose Fraction’) 
ylabel(‘Mid-body Fraction’) 
Zlabel(Yvprime Percent Error’) 
axis([.05 .25 .40 .60 -20 20]) 
view(60,30) 
print -depsc2 yvperr 
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